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The breathing-mode giant monopole resonance is studied within the framework of the rel- 
ativistic mean-field (RMF) theory. Using a broad range of parameter sets, a systematic 



analysis of constrained incompressibility and excitation energy of isoscalar monopole states 
in finite nuclei is performed. A comparison is made with the incompressibility derived from 

ON ■ 

the semi-infinite nuclear matter and with constraint nonrelativistic Skyrme Hartree-Fock cal- 
culations. Investigating the dependence of the breathing-mode energy on the nuclear matter 

1 : tocom pressibility ' " 18 8how " that dynamical properties of surface respond d,fferently in the 

■ RMF theory than in the Skyrme approach. 
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The breathing-mode giant monopole resonace (GMR) has been a matter of contention 
in the recent past ffl. The energy of the GMR has been considered to be a source of 
information on the nuclear matter compressibility. Empirically, data on the GMR have 
been measured with considerable precision 0. The detailed analysis of the GMR leading 
to the incompressibility of nuclear matter is as yet not settled and is still under investiga- 
tion. Theoretically, the incompressibility is understood to have been obtained || using the 
density-dependent Skyrme interactions. The deductions base themselves upon interpolation 
between various Skyrme forces for the GMR energies obtained from HF+RPA calculations. 
There has been, however, no unambiguous reproduction of the empirical GMR energies of 
medium-heavy nuclei using this approach. The breathing-mode energies depend not only 
upon the bulk incompressibility, but are also sensitive to the surface incompressibility. The 
relationship between the bulk and the surface incompressibility has been enunciated clearly 
for the Skyrme approach in Ref. ||. 

The relativistic mean-field (RMF) theory has achieved a considerable success || in de- 
scribing the ground-state properties of nuclei at and far away from the stability line. The 
dynamical aspects have, however, remained largely unexplored. The breathing-mode ener- 
gies and incompressibilities were obtained within the RMF theory using the linear Walecka 
model PJ^]. The relationship of the GMR energies to the incompressibility of nuclear matter 
is not yet clear for the RMF theory. In the non-relativistic Skyrme approach, on the other 
hand, the GMR energies are connected to the incompressibility straightforwardly. In this 
letter, we examine the dependence of the breathing-mode energies in finite nuclei on the 
incompressibility of nuclear matter in the RMF theory. 

We start from a relativistic Lagrangian || which treates nucleons as Dirac spinors ip 
interacting by the exchange of several mesons: scalar a-mesons that produce a strong at- 
traction, isoscalar vector cu-mesons that cause a strong repulsion and isovector p-mesons 
required to describe the isospin asymmetry. Photons provide the necessary electromagnetic 
interaction. The model Lagrangian density is: 

£ = ^{ij^ - M}i/> + \d^ad^a - U(a) - g^cnp 

where U(a) is the non- linear potential with the cubic and quartic terms [BJ: 

U{a) = ^m a a 2 + ^cr 3 + ^g 3 a 4 . (2) 

M, 77V, vtiu and m p are the nucleon, the a-, the u-, and the p-meson masses, respectively, 
and g a , g^, g p and e 2 /47r=l/137 are the coupling constants for the a-, the u-, the p-mesons 
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and for the photon. The field tensors for the vector mesons are Vt^ v = d^u v — d u u p , R^ u = 
d^p u - d v p p - g p (p» x p u ) and for the electromagnetic field = d»A v - d u A p . The 
variational principle leads to the stationary Dirac equation with the single-particle energies 
as eigenvalues, 

h D ipi(x) = ei^i(x), (3) 

where 

h D = -id. V +P(M + g a <j(r)) + g u u°(r) + g P r 3 p°(r) + e ^A°(r). (4) 

Solving these equations self-consistently one obtains the nuclear ground state in terms of 
the solution 

In order to obtain the isoscalar monopole states in nuclei we perform constrained RMF 
calculations solving the Dirac equation 

(JiD — Ar 2 ) ipi(x, A) = eiipi(x, A), (5) 

for different values of the Lagrange multiplier A which keeps the nuclear rms radius fixed at 
its particular value 

R={\jr"p x (r)d 3 r] 1/2 , (6) 

where 

A 

px(r) = p v (r, A) = J2 ^l( x > Wi( x > A) (7) 

i=l 

is the local baryon density determined by the solution ipi(x, A). The total energy of the 
constrained system 

E RMF (X)=E RMF yj i (x,X)}, (8) 

is a function of A (or the rms radius R) which has a minimum, the ground state energy 
Ermf = Ermf{Q), at A = corresponding to the ground-state rms radius Rq. 

This behaviour of the constrained energy (8) as a function of A (or R) allows us to 
examine the isoscalar monopole motion of a nucleus as harmonic (breathing) vibrations 
changing the rms radius R around its ground-state value Rq- Considering s = (R/Rq — 1) 
as a dynamical collective variable and expanding (8) around the ground-state point s=0 (or 
A = 0) we obtain in the harmonic approximation, 

E RMF (X) = E° RMF + l -AK c (A)s\ (9) 
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where 



K C {A) = A' 1 [ R 



,d 2 E 



RMF 

dR 2 



(A)' 



(10) 



A=0 



is the constrained incompressibility (TD| of the finite nucleus. The second term in eq.(9) 
represents the restoring force of the monopole vibration. In order to obtain its associated 
inertial parameter we apply the method |J of making a local Lorentz boost on the con- 
strained spinors ipi(x, A): 



ipi(x,t) = —S(v)ipi(x,X), 
where S(v) is the Hermitian local Lorentz boost operator 



(11) 



S(v) = cosh ( — 1 + a.pr sinh ( — 



(12) 



with cf) = tanh 1 (\v\) and 7 = Vl — v 2 - The velocity field v = v(f,t) is then obtained by 
the continuity equation 



d 



dt j 



Yjl>\( x , ^)M x i A) + V-H^( x ' A)a?M£, A) = 0, 



(13) 



i=i 



which, using eqs.(7) and (11), transforms into the form 

. dp(r,X) 



ds 



+ v.(p(r,\)v(r,t)) = 0. 



(14) 



Up to first order in s the velocity field is determined by eq.(14) into the form v = —su{r)r/\f\ 
where 



u(r) 



'0 as 
p(r, A)r 2 



(15) 



The inertial parameter of the monopole vibration is then obtained as 

B rel (A) = A- 1 J u 2 (r) H RMF (r) d 3 r, 



(16) 



where u(r) is the velocity function (15) at A = and TCrmf(^) is the Hamiltonian density. 

We have obtained K C (A) from eq.(10), B rd {A) from eq.(16) and the frequency of the 
isoscalar monopole vibration as 



K C {A) 

\ Brel(A) 



(17) 
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for a number of spherical nuclei. Various parameter sets such as NL1 ||11|| , NL-SH ||, NL2 



1~2|| , HS |T3| and LI [[HJ with values of the nuclear matter incompressibility K^m =211.7, 



354.95, 399.2, 545 and 626.3 MeV, respectively, have been employed in the calculations. 
The last two sets, HS and LI, correspond to linear models without the self-coupling of the 
cr-field. In addition, the set LI excludes the contribution from the p-field. Among the sets 
NL1, NL-SH and NL2 which correspond to the non-linear model, only the set NL2 has 
a positive coupling constant g 3 , eq(2). Whereas the set NL1 reproduces the ground-state 
properties of nuclei only close to the stability line due to the very large asymmetry energy, 



the set NL-SH describes also nuclei very far away from the stability line |L4 . 

Results from the present constrained RMF calculations are shown in Figs. 1 and 2, where 
the collective mass B re i(A), the constrained incompressibility Kq(A) and the associated 
excitation energies uq for a few nuclei have been displayed. First, we consider Kq{A) in 
Fig. 1 (a). The incompressibility of nuclei shows a strong dependence on the nuclear matter 
incompressibility K^m, with a few exception for light nuclei. For the linear force HS, Kq 
shows a slight dip from the increasing trend for 208 Pb and 90 Zr, whereas for light nuclei 
40 Ca and 16 the HS values are even smaller than the NL2 values. The dependence of 
the imcompressibility K(A) of finite nuclei on K^m obtained from non-relativistic Skyrme 
calculations is different: it increases monotonically with K^m 0. This difference can be 
understood from the difference in the behaviour of the surface incompressibility in the two 
methods. In the Skyrme approach, the surface incompressibility has been shown to be 
Ks ~ —Knm for standard Skyrme forces ||. However, this does not seem to be the case for 
the RMF theory as shown by the dip at the HS values. Thus, the surface incompressibility 
is not necessarily a straight function of the nuclear matter incompressibilty in the RMF 
theory. 

Since the GMR energy depends strongly upon the inertial mass parameter, we next 
examine the mass parameter B re i(A) in Fig. 1 (b). The collective mass obtained in Eq. 16 
from the velocity field shows an interesting dependence on the mass of the nucleus considered. 
For the heavy nucleus 208 Pb, the collective mass is about 0.9 and for 90 Zr it is 0.5 MeV" 1 . For 
lighter nuclei it decreases to about 0.3 MeV -1 , as shown in Table 1. It does not, however, 
depend much upon the parameter set used and therefore shows only little sensitivity to 
Knm- This implies that the Knm dependence of the energy is then predominantly due to 
the K NM dependence of the incompressibility K C (A). 

We now compare in Table 1 the results from two RMF parameter sets, NL1 and NL-SH. 
The constrained incompressibilities K~c for NL-SH are higher than for NL1 as discussed 
above in Fig. 1(a). It is interesting to compare the relativistic collective mass B re i(A) with 
the expression usually applied in the non-relativistic sum-rule approach, 

B sr {A) = MRl (18) 

where the ground-state rms radius Rq is used. The ratio B sr (A)/B re i(A) is shown in Table 
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1 where the energy uc, eq.(17) and its approximation 



have also been compared for the sets NL1 and NL-SH. The observation that the collective 
mass B re i(A) decreases significantly (up to 50 %) for light nuclei, provides a hint for the 
influence of the surface in the value B re i(A). For heavy nuclei the approximate value B sr (A) 
is rather close to B re i(A). Consequently, the monopole energies ujc and u sr differ by about 
1 MeV only for heavy nuclei and up to about 5 MeV for the lighter ones. 

The energy uj sr corresponds to the monopole excitation energy E\ usually obtained from 
nonrelativistic Skyrme forces within the sum rule approach |Tj| and in nonrelativistic con- 
strained Hartree-Fock calculations. In Table 1 energies ujq and u sr have been compared with 
such nonrelativistic constrained HF results tosky obtained from Eq. (19) with the Skyrme- 
type forces SkM and SIII having about the same nuclear matter in compressibility K^m as 
the sets NL1 and NL-SH, respectively. From Table 1 it can be seen that the Skyrme results 
J|] u>sky are actually close to the energies u sr and differ significantly from the values of u c . 
This difference in the RMF constrained energy ujc from the Skyrme constrained energy uJsky 
is small for heavy nuclei. It, however, increases for lighter nuclei, where the RMF shows 
smaller values. The difference is particularly significant for the higher compressibility forces 
(NL-SH and SIII). It arises naturally from the lower values of the collective mass B re i(A) 
for light nuclei in the RMF theory as discussed above. Thus, the masses B re i and B sr . have 
different dependences on surface. 

Fig. 2 shows the constrained breathing-mode energy ujc for different nuclei. For heavier 
nuclei 208 Pb and 90 Zr, ujc shows a behaviour similar to that shown by Kc in Fig. 1 (a), 
as the collective mass does not show much change with incompressibility (Fig. l.b). The 
stagnation in energy at NL2 and HS is reminiscent of the effect of surface compression as in 
Fig. 1 (a), thus reflecting the role played by the surface in the RMF theory. In lighter nuclei 
this effect is more transparent. For 40 Ca and 16 0, ujq is not related to K^m i n a simple way 
due to the combined effect of Kc and the collective mass B re i, where there is even a decrease 
in u>c with K NM . Employing schematic parameter sets, it was also shown [|T?|] earlier that 
within the relativistic quantum hadrodynamics the surface compression responds differently 
than in the Skyrme ansatz. 

A comparison of the empirical values of the GMR energies with uoc is worthwhile. The 
empirical values for 208 Pb and 90 Zr are 13.9±0.3 and 16.4±0.4 MeV, respectively. For lighter 
nuclei the values are very uncertain. Systematics of the values for 208 Pb show that the 
empirical value is encompassed by the constrained calculations curve from K n m = 300 - 400 
MeV. The 90 Zr value is, however, overestimated by the corresponding curve by 2-3 MeV. 

We have also carried out calculations of the incompressibility for finite nuclei by scaling 
the density using semi-infinite nuclear matter. The incompressibility for a finite nucleus 
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can be written as: 

K A = K NM + KsA- 1 " + (^-) 2 + | - [ 1 - 27^-] ■ ( 20 ) 

\ A J 5r \ K NM J 

where major contribution to Ka arises from the volume {K^m) and the surface terms (Ks). 
The surface incompressibility Kg has been obtained by calculating the second derivative of 
the surface tension (a) with respect to the density for each change in the scaled density as 
given by @ 

K s = 4vrr 2 {22a( Po ) + 9^"(po) + ^^P^(Po)}, (21) 



The details of the procedure to perform scaling of the density have been discussed in ref |L6| . 
The last two terms in eq. (20) contribute very little. We add these terms for completeness, 
however. The asymmetry coefficient K^ has been taken at —300 MeV from the empirical 



determination [18] and is a reasonable value. The third derivative of the EOS, p^e'", for each 
force is known from the nuclear matter calculations. We have performed the calculations 
with the parameter sets used above. 

Fig. 3 shows the 'scaling' incompressibility Ka for various RMF forces obtained from 
Eq. (20). The general trend of the scaling incompressibility with Knm is about the same 
as in Fig. 1(a), showing a dip for the force HS. In general, the Ka values are about 20% 
larger than the constrained values Kc- This is consistent with the known relationship of 
the two types of the incompressibilities for nuclear matter, whereby it was shown [|10j that 



Ka(NM) = ^K C (NM). Fig. 3 brings out again the importance of the role of the surface in 
the dynamical calculations within the RMF theory. The surface incompressibility Ks for the 
forces NL1 and NL-SH are —333.1 and —610.1 MeV respectively. Thus, in both the cases the 
ratio of the surface to the bulk incompressibility is obtained as 1.58 and 1.72 respectively. 
These values differ considerably from the ratio of 1 in the Skyrme ansatz. The values of the 
bulk and surface incompressibilities for NL-SH are closer to the empirical values from Ref. 

§■ 

In conclusion, it has been shown that within the RMF theory, the incompressibility 
of a finite nucleus, K(A), does not depend on the nuclear matter incompressibility in a 
simple way. This dependence is strongly related to the properties of the surface subjected to 
compression. This is different from the behaviour of the surface in the Skyrme ansatz. The 
collective mass for the breathing mode vibration in the RMF theory also shows a different 
behavior from light to heavy nuclei, thus affecting the frequency of the isoscalar breathing 
mode. 
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TABLES 

ABLE I. The constrained incompressibility Kc(A) in MeV, the mass parameter B re i(A) in 
, the ratio B sr (A) / B re i(A) and the associated monopole frequencies u>c and ui sr , (both in 
) calculated with the sets NL1 and NL-SH. Comparison is made with constrained Skyrme 
.ts oJsky obtained within the sum-rule approach |l5| with the Skyrme forces SkM and SIII. 

NL1: K NM = 211.7 MeV NL-SH: K NM = 354.95 MeV 

SkM: K NM = 216.7 MeV SIII: K NM = 356.00 MeV 



ei K c(A) B(A) B sr jB Td lo c io sr oj S k y K C {A) B(A) B sr /B re i w c uj sr Ld S k y 
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FIGURES 

'IG. 1. (a) The constrained incompressibility Kc in Eq. (10) for a few nuclei obtained using 
•us RMF parameter sets, (b) The collective mass B re i for the breathing-mode monopole 
itions obtained from Eq. (16) within the RMF theory. 

FIG. 2. The frequency u>c of the monopole mode obtained using eq. (17). 

'IG. 3. The incompressibility K\ (Eq. 20) obtained from 'scaling' of the nuclear density in 
;emi-infinite nuclear matter using the Thomas-Fermi approximation. 
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